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ABSTRACT 



Context. The mass function of super-massive black holes in our cosmic neighborhood is required to understand the statistics of their 
activity and consequently the origin of the ultra high energy particles. 

Aims. We determine a mass function of black hole candidates from the entire sky except for the Galactic plane. 
Mettiods. Using the 2MASS catalogue as a starting point, and the well established correlation between black hole mass and the bulge 
of old population of stars, we derive a list of nearby black hole candidates within the redshift range z < 0.025, then do a further 
selection based on the Hubble-type, and give this as a catalogue elsewhere. The final list of black hole candidates above a mass of 
Mbh > 3 • 10' Mq has 5,829 entries; moreover doing a further Hubble type correction to account for the selection effects cuts down 
the number to 2,919 black hole candidates. Here we use this catalogue to derive the black hole mass function. We also correct for 
volume, so that this mass function is a volume limited distribution to redshift 0.025. 

Results. The differential mass function of nearby black hole candidates is a curved function, with a straight simple power-law of index 
-3 above 10* Mq, growing progressively flatter towards lower masses, turning off towards a gap below 3 10* Mq, and then extending 
into the range where nuclear star clusters replace black holes. The shape of this mass function can be explained in a simple merger 
picture. Integrating this mass function over the redshift range, from which it has been derived, gives a total number of black holes 
with z < 0.025, and Mbh > 10' M© of about 2.4 • 10*, or, if we just average uniformly, 0.6 for every square degree on the sky. 
In some models many of these, if not all, are candidates for ultra high energy particles sources. Even if a very small fraction of the 
super-massive black holes produces ultra high energy cosmic rays, this should be enough to observe from the statistics of their arrivals 
the highly inhomogeneous distribution of the galaxies and their super-massive black holes; also this distribution may be smeared out 
by the possible scattering of the arrival directions of the particles by intergalactic magnetic fields. 

Key words. Black hole physics - Galaxies: general - Acceleration of particles 



1. Introduction 

In the quest to find the origin of ultra high energy cosmic ray particles, many proposed mechanisms use the power of super-massive 
black holes, either obtained from accretion, e.g. ( Falcke & Biermann|1995 1, or from spin-down (e.g. ( Blandford & Znajek|1977[ l), 



or in some other way (e.g. (Biermann & Fram pton|2006| l). It is not a priori clear, which is the best constraining argument one should 
use to find out what possibly very small fraction of all super-massive black holes do produce ultra high energy particles. As long 
as the mechanism is not clear, we need to establish a basis for discussing all plausible super-massive black holes. Therefore it is 
important to obtain the mass function of black holes in order to estimate the number of black holes to within the sphere from where 
ultra high energy particles might come. 

Since various other groups have also derived the mass function of black holes, using different methods, we have a good way to 
check our work. 

Super-massive black holes characteristics are strongly correlated with the propeities of the host galaxies, and we use this coiTe- 
lation to obtain abundant statistics, but also to control the quality of our data set. We wish to address a number of questions related to 
the super-massive black holes in the universe: a) It appears that super-massive black holes have a low cut-off in the mass distribution, 
and it is important to verify this and try to understand why there is such a characteristic mass; also, it appears that the black holes' 
mass function has an upper limit, and again we need to understand whether this is just a statistical fluke, or has a physical meaning, 
b) What is the slope of the mass function? This will surely constrain our idea, on how black holes are growing, c) What is the energy 
input into the universe during the growth? This is an additional constraint, d) What is the number of black holes? This will constrain 
any search for directional correlations with ultra high energy cosmic ray events. 
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2. Derivation 

Super-massive black holes are common in the centers of galaxies, and their mass is correlated with the properties of the surrounding 
galaxy, for instance with the velocity dispersion of the nearby stars, and with the mass of stars in the spheroidal distribution of the 
old stellar population. 

Since the old stellar population has a spectrum peaking near 2 /i {vSy) (van der Wei et al. 2006), and it is this population, that 
correlates with black holes, we use the 2MASS catalogue as our starting point ( |Skrutskie et al.,2006): We focus on e arly Hubble 
type galaxi es (jCaramete & Bierm ann 20 10|, and us e the black hole mass spheroid correlation, e.g. ( Faber et al.|1997| Magorriali| 



etaLl|T998l |Wang & B iermann l998f |Silk & R ees 1998] [Haring & Rixl|2004{ [Tremaine et aL][2052 



Aller & Richstone||2007| ) to 



derive the black hole mass function. We select down to 0.03 Jy at 2 This is far above 2MASS catalogue completeness limit. 
The procedure is as follows: 

Step 1: We have first the 2 micron flux density limited sample of 10,284. We use all the distance corrections available for this 
sample. 

Step 2: We limited the sample to well determined Hubble types, and only early Hubble types, and have 5,894. 
Step 3: The decisive step is to eliminate all galaxies, for which we did not have sensitivity for the presence of a black hole of a 
given mass equally for all Hubble types. This is a decisive step, and cuts down the sample by another factor of two, to 2,928. 

For this sample of 2,928 we show the luminosity function that shows the classical Press-Schechter function ( Schechter|1976| l, 



a power-law with an exponential cutoff. This agrees with what is available in the literature, but uses a different approach, and for 
black holes very much larger numbers. 

Step 4; We get the semi-final BH candidate mass function, extending apparently down to about 10^ solar masses. 
In detail then: 

The mass of the central black hole is proportional to the flux density of the stellar emission at 2 /i and also propor- 
tional to the luminosity distance squared. We need to determine the proportionality constant. This constant is, however, dif- 
ferent for different Hubble types of galaxies, and we have to determine so the constant separately, and then check for self- 
consistency. Since the spheroidal stellar population is much smaller in later Hubble types, these differences can be large. Because 
the estimate of the black hole masses uses the redshift, which is a bad indicator of true distance at distances less than that 
of the Virgo cluster, we use the available distances provided by the work of Barry F. Madore and Ian P. Steer (available at 
http://nedwww.ipac.caltech.edu/level5/NEDlD/intro.html) who compiled a database of 3,065 accurate, contemporary distances to 
1,073 galaxies with modest recessional velocities (that is, less than 1/8 c) published between almost exclusively 1990 and 2006. 
From this we match 429 distances in the catalog of massive black hole, and for the rest of distances less than that of the Virgo cluster 
we used distances with respect to the Virgo infall only provided by NED. 

We use known black hole masses to calibrate the various early Hubble types; this control sample comprises 58 black holes 
selected from the recent literature; we use only galaxies earlier than type Sc, including type Sb. These fits give us an error bar of 
± .4 dex, slightly worse tha n the ± 0.3 dex error given by more sophisticated fits using a much better but also much smaller data 
set ( Aller & Richstone|2007 ). These corrections affect a large fraction of the candidate black holes. The final list above 3 ■ 10^ M© 



comprises 5,829 black hole candidates. The main error is + 0.4 dex, and affects most of all the binning of the sample. 

In table [T]we show, how many galaxies of various kinds are excluded, and it becomes obvious, that from these statistics, at 
worst we have a factor of less than 2 error, if all excluded galaxies of no known Hubble type turned out to be E or SO. 



Table 1. 



Type 


5894 selection 


2928 selection 


Rejected list 


E 


783 


765 




SO 


2626 


1771 




Sa 


879 


184 




Sb 


1052 


128 




Sab 


554 


80 




Sc 






1941 


h-r 






14 


Sy 






128 


No type 






867 


Sbrst 






421 


Unknown type 






1019 



We show the correction factor dependent of the Hubble type, relative to elliptical galaxies in Fig. [T] 

In histogram |2]we show that we have galaxies with black hole masses over the entire range of to 10^ as calibrators. 

We derive a mass function by going to the minimum mass at each redshift, for which the list is complete to the chosen flux 
density in the 2MASS catalogue, starting with 0.025 and going down, always correcting for smaller volume. The sequence of steps 
is as follows: we obtain the number of black holes in several bins above the minimum mass we can obtain for both spirals and 
ellipticals, at redshift z - 0.025; then we go down in redshift, correcting for different accessible volume for a complete sample. 
Figure |3]illustrates the procedure. 

This diagram of black hole mass versus redshift shows clearly, that near the flux density limit, at any redshift, we have only 
spirals, since for a given redshift and given flux density, elliptical galaxies have an inferred black hole mass a factor of up to 5 larger 
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Fig. 1. This is the correction factor to obtain the black hole mass from the 2 yU emission of the old stellar population, as a function of 
Hubble type, relative to elliptical galaxies. 
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Fig. 2. Histogram of the calculated black hole masses for the sample of black holes obtained from the Uterature. 

than spirals, Sb spirals being the most extreme. In order to eliminate this bias, we have cut the entire sample again such that we 
have a black hole mass limit at each redshift, which allows both spirals and ellipticals. This procedure eliminates many, but not all 
of the late type spirals from our final sample to determine the mass function. The structures seen in this diagram illustrate the bias 
introduced by the correction factor to derive black hole masses from flux density limited data. This cuts down the fraction of spirals 
in the final "Hubble-selected" sample to about twenty percent, so reducing greatly the error contribution from the Hubble-type 
correction shown above. 

Therefore we follow the curved red line in Fig. [3] get the number in that mass bin, and ratio to volume; and so we keep adding 
bins at lower mass until we reach the smallest bin, we wish to consider. This procedure ends at redshift z = 0.0045 for M > 10^ Mq, 
and at z = 0.0024 for M > 3 ■ lO*" Mg. 

Next we show the luminosity and mass functions, using this final selection of galaxies and their black hole candidates. 

The luminosity functions (Fig. |4| correspond to what is in the literature, as we discuss further below; they are all well fitted 
with a Press-Schechter law (1974), a power-law with an exponential cutoff. Since our completeness criterion is different from what 
most others have used, we note that our luminosity functions are similar Since we selected far above the completeness limit of the 
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Fig. 3. Plot of Mass over redshift for the massive black hole catalog with the selection curve for elliptical galaxies in blue. This also 
shows the stepwise selection procedure with the two most extreme cases, the lowest mass, and the highest redshift, also in red. 
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Fig. 4. Integral luminosity function corrected for Hubble type sampling, 2928 objects selected in K-band, the parameters of the 
Schechter function fit are: M. = -24.3, 0. = 0.0008, a = -1.8. 



2MASS sky survey, and aim to almost 80 percent at elliptical and SO-galaxies, our procedure of selection may be more complete 
for the task at hand, to determine the black hole mass function, independent of any other property of the black hole or the activity 
of its surroundings. 

The errors included here include an estimate from the binning error induced by uncertain distances, and of course, also Poisson 
statistics; these are really relevant only for the highest mass bin, where we have 25 galaxies with their presumed black holes. 
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increases our 



We have consciously ignored all super-massive black holes in late Hubble type galaxies; this underestimates our numbers. 
Correcting for the Galactic plane using a histogram of numbers per latitude bin (see Caramete & Biermann||2010 1 
number by 12.5 percent. 

We show the result in Fig. |5] with error bars and fits through the data. 
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Fig. 5. Integral mass function corrected for Hubble type sampling, 2928 objects, the slope of the lines is: thick line -2.0 fitting 
> 10*^ Mq, and thin line -1.0 fitting between 10^ and 10^ Mq. 



Since it is unlikely that all these galaxies actually include a central black hole, we cut the distribution near 3-10* Mq, the mass 
of our Galactic Center black hole; we obtain, including the correction for the missing Galactic plane (12.5 percent) 



(p(> Mbh) = 6-10 



-3±().4 / MBH,mm 



1.0 



Mpc 



between about 10^ Mq and about 10^ Mq, and 



0(> Mbh) = 9-10- 



4±0.4 / ^BH.min 

10«Mo 



-2.0 



Mpc-\ 



(1) 



(2) 



above about 10^ Mq. These two laws are reasonably good fits, and match near 1.5 ■ 10^ Mq. We imposed the exponents pragmat- 
ically, since they can be determined from the data only with large error bars, strongly depending on which bins are included. Simple 
power-laws are adequate. 

Between a nominal 10^ Mq and 3 ■ lO*" Mq we have only 6 candidates, also between 3 ■ 10^ Mq and 10^ Mq we have 46 
candidates, and between 10^ Mq and 3 ■ 10^ Mq we already have 185. Also, in the mass range between 10* Mq and 10^ Mq there is 
increasing doubt with lower masses, whether we really always have a black hole, and not just a nuclear star cluster. 

Comparing this and various other fits which we have tried, we summarize these attempts including error bars from the original 
calibration to various Hubble types; however, after our correction for Hubble type selection there are few Sb spirals left in the sample 
(128). Comparing a fit through only the ellipticals and SO's and the entire sample, with a cut of low mass, and without gives us an 
estimate of the error in the absolute normalization of about 0.25 in the log; since there are systematics in allocating Hubble types, 
and including or excluding galaxies, we increase our estimate of the error by about V3 to get 0.4 in the log. We estimated the errors 
on the absolute scaling, and the exponent from the fit. A lowering of the low mass cutoff of a factor of 3 from lO'^ Mq to 3 ■ lO*" Mq 
corresponds to a negligible increase in the integral number; the integral mass function is nearly flat in that range, which of course 
implies that the differential mass function in this interval is consistent with zero. A final test was to cut in the mass-redshift diagram 
even more severely to check on edge effects in the sampling; the black hole density increased marginally, less than 0.2 in the log. 
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As any evolution of the black hole mass function is expected to be approximately with (1 + z)'^'^^, up to a redshift of 0.025 any 
evolutionary effects ought to be negligible. On the other hand, it is not clear, that we have reached the distance range, over which 
the universe is already perfectly homogeneous: that distance might be larger than 500 Mpc (Rudnick et al. 2007, Kashlinsky et al. 
2008); there is no evidence that the galaxies in the universe change their statistics drastically between 100 Mpc and 500 Mpc or 
even more of the local universe. 



300 



250 




200 



^ 150 



100 



50 




0.005 



0.010 



0.015 



0.020 



0.025 



Fig. 6. Histograms of the black hole candidates, corrected for Hubble type sampling, Hubble type correction, cut in BH mass at 
1 * IO^Mq, 2023 objects, with a fit valid for a homogeneous universe. 

In Fig. |6]we show our numbers of black holes above 10^ Mq as a function of redshift and the fit up to the redshift of a little above 
0.015, where our procedure should be fairly complete for such a mass limit. A homogeneous universe fit is amazingly adequate, but 
the data also clearly demonstrate that near the redshift limit we are beginning to miss galaxies, here about 15 percent of the total. 

3. Comparisons 

There are a large number of explorations of the black hole mass functions, and their energy input, and here we compare with those 
papers. 



3. 1 . Numbers of black holes 

Obviously, we cannot be certain tha t actually all these black hole candidates are really black holes: this question was thoroughly 
investigated by [Gree ne et al. (j2008j), who could not find a large number of black holes below 10^ Mq. What we cannot easily say, 
is to what degree this amounts to a selection effect. On the other hand, it has been shown by Cote et aL] p006 ), that the correlations 
between galaxy parameters and the super-massive black hole apparently seamlessly merge into a correlation with nuclear star 
clusters. Nuclear star clusters can in turn become super-massive stars, which collapse to black holes ( |Sanders|1970[[Portegies Zwart| 
et al.|200 4 ). This may in turn also explain the mass of the black holes near the transition. 

The mass function is probably best determined near a ma ss of 10 *^ Mo, where we obtain 9 ■ lO""*"^"'*" Mpc^^. At that mass the 
density is in agree ment, within the error bars, with bothlLauer et al. 62007), who show a possible range from 2 ■ 10^'^ Mpc^^ to 



5 ■ IQ-^ Mpc-\ and 



Benson et al. 



( |2007[ ), who show a number near of 2 ■ lO^^* Mpc^^ (correcting for Hubble constant). 
Greene & Ho ( 2007a| l have determined the black hole mass function, emphasizing the low masses, and selecting for activity: 
First, their numbers show, that this distribution peaks near 10^ Mq, confirming that it may be reasonable to cut at that mass. Second, 
they confirm the slope which we obtain. However, their total numbers are very much lower - by a factor of about 5000 - than in 



Lauer et al. ( 2007 1; Benson et al. 



based on the activity of black ho 



( 2007 1, and our results, stemming from their very different selection procedure; their selection was 



es. 



L.I. Caramete and P.L. Biermann: The mass function of nearby black hole candidates 



7 



We conclude here, that our statistics are consistent with eariier numbers. 



3.2. Energy input from black holes 



As black holes grow, they accrete from an emitting accretion disk; the luminosity depends on the growth rate, and so the final mass 
is a direct measure of all growth (Liist 1952; Sh akura & Sunyaev| |1973 ; Noviko v" & Thorne| [l973 ). Similarly, all well measured 
active galactic nuclei have radio emission, e.g., ( |Nagar et al.|2 005C which we interpret to come from a relativistic jet, e.g., ( Falcke| 
|& Biermann||1995j ); however, also in this case, all energy dissipated and emitted by the jet is a measure of growth. On the other 
hand, mergers of black holes, e.g. ( Gergely & Biermann|2009 1, would more likely emit abundant gravitational waves, invisible so 



far in the normal universe. In order to study this we need to obtain the total growth in black holes, so their total mass. 
If we integrate from 10^ Mq, the total mass density in black holes is given by 

2.8 ■ 10^^° '* Mo Mpc-^ (3) 

barely dependent on where we cut the distribution, since at low masses the integral function goes flat, and at high masses the 
distribution goes steep. 



This compares, e.g., with the number obtained by Kronberg et al. ( 2001 1, which has been derived from radio emission signatures, 



and so depends on the efficiency of the conversion process: they obtain a lower limit of 

; 2.2 ■ 10^ e(xJ Mo Mpc-\ (4) 

where eo.i is the efficiency in units of ten percent. [Kronberg et al.| ( [200T] ) derived their numbers from a selection of higher mass 
black holes, and give a lower limit. One could argue that all power output by black holes in the mass range between 10^ Mo and 
10^ Mo influences mostly the host galaxy, and is less relevant for the intergalactic medium; it could, for instance, mostly contribute 
to driving a galactic wind (see [Everett et al . 2008; Breitschwerdt 2008). Above we ignored this contribution. 

Another way of looking at this is by simply converting the units, and then the mass density derived from this is 

1.9-10-™°gcm-3 (5) 

which corresponds to a cosmological density of (using for the Hubble constant 70. 1 km/s/Mpc, the critical density is 9.77 ■ 10"^*^; 
[Komatsuet al.X2009j 

^BH = 2- l0-«*°-40 (6) 
where the error bar in the first factor corresponds to a range of 0.8 ■ 10"^ to 5 ■ 10"*. 

Assuming, that growth of black holes is mostly via baryonic accretion ( [Shakura & Sunyaev|1973 Novikov & Thorne 1973 1, the 



energy input into intergalactic space from the production of all these black holes is given by 
^ 9. io-i5±0-40^yjg^^„„-3^ (7) 

which is by a factor of order 4.5 ■ 10""^ - excluding the uncertainties - far below that due to the microwave background (at 
4.2 ■ 10"'^ erg cm^^); this number is an average over all space, derived from the local universe to redshift 0.025. 

We need to compare this estimate with the energy density visible directly in electromagnetic background, with the energy density 
inferred to be there in the form of turbulent motion, heat, magnetic fields, energetic particles, as well as comprehensive simulations 
of all such processes. 

The electromagnetic radiation background in other wavelengths is approximately (compilation by Kneiske (2009), private com- 
munication): at FIR and optical frequencies lower than the microwave background by 5 ■ lO^^*, in X-rays down by 5 ■ 10"^, decreasing 
with higher frequencies then to a factor of 1.5 ■ 10 * down for y-rays. The energy input from the growth of black holes, even at only 
10 percent efficiency as used above, exceeds all other energy densities, and could be marginally consistent with the FIR and optical 
background. 

The energy input has also been estimated by [Gopal-Krishna & Wiita| ( [2001 [ ), and [EnBUn et al.[ ( [1998| l. EnBlin et al. obtain 
(2.5 -4.0) ■ 10" ""erg/cm^, while Gopal-Krishna et al. get 2 ■ IQ^^^erg/cm^ in filaments, similar to what we obtain. We conclude that 
our numbers are consistent with earlier estimates, derived differently. 

We will compare with the very uncertain energy densities of ultra high energy particles elsewhere, since such a comparison 
strongly depends on the spectrum of energetic particles below the energies at which we can observe them. 

We should also compare this with the numbers derived from simulating the turbulent dynamo in the intergalactic medium ( |Ryu 



|et al.|[2008] l; Ryu & Kang (2009), private communication: in this case we need to compare with the energy input into thermal 
motions and magnetic fields: They obtain a magnetic field energy density of 2 ■ 10"'^ erg/cm^ averaged over the universe at zero 
redshift; in filaments a magnetic field energy density of 1 ■ 10"'* erg/cm^, in average a thermal energy density 3 ■ 10"'* erg/cm^, and 
in kinetic energy a density 1 ■ 10"'^ erg/cm^. All this energy in these simulations derives from the gravitational wells of the dark 
matter large scale structure, including the galactic halos. To inject far more energy from the activity and growth of central black 
holes could present a serious problem. The energy input from active galactic nuclei into magnetic fields might be far less than 10 
percent efficient; to satisfy observational constraints we have the condition then 

1 ■ 10-2±0-4 > (8) 
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where eg is here the efficiency to channel energy from black hole growth into overall magnetic field. We might satisfy this by 
arguing that much less than 10 percent of all black holes ever become radio galaxies, or alternatively that all super-massive black 
holes are radio galaxies for much less than 10 percent of their time, and much less than 1 percent of their accretion energy is 
charmelled into magnetic fields (we need a factor of 2 ■ 10"^). This may then be consistent with the reasoning of |Kronberg et al.| 
( 2001| l, who use just the manifest energy blown into the intergalactic medium by radio galaxies, traceable through their nonthermal 
radio emission, to derive numbers for the black hole density. 

There are several ways out of this possible dilemma: an obvious one is that the fraction of energy going into various channels, 
electromagnetic, kinetic, magnetic fields is far from settled. A last resort is to note that much of the energy emitted during the growth 
of black holes via mergers could go into gravitational waves. 



3.3. Extrapolation to lower masses? 

Now, we have to remember that we cannot be absolutely sure that the low mass cutoff is in fact at 10^ Mq, it might be at stellar scale, 
e.g., ( [Clark et al.|2002[ Irwin et al. 2009 ). The re are very few confirmed black hole candidates with a mass, that significantly exceeds 



10 Mq, and is below 10^ Mq: (see|van Kerkwijk et al. 1996; Trudolyubov & Priedhorsky 



2004; Gebhardt et 



al.|2005 



Patruno et al. 



2006} [Maccarone & Koerding||2066[ JOrosz et al.||2007, ,Ulvestad et al.||2007[ [Casella et al.,,2008, ,Siiverman & Filippenko||2008 



Abubekerov et al.|2009t|Farreli et al.|2009| l. Oneof the most convincing cases for an intermediate mass black hole is the candidate 
in the cluster Gl, galaxy M31 also known as Andromeda. 

Since the integral mass function flattens out just near the mass, where there may be very few black holes anymore, below 
3 ■ 10^ Mq, there maybe little mass there. The mass to there is only 2.8 lO^'^" '* Mq Mpc^^. This is far below the numbers for stellar 
mass black holes, according to some estimates ( |Agol & Kamionkowski|2002| ; it is intriguing to note that the total mass in stellar 
size black holes inside our Galaxy is estimated to be of order a few 10 Mq- '^'"^ ^^^^rW '^r^'' ^rvnnor-t trv this it ic tVion 
reassuring. 



2006 



, our law clearly does not connect to this. It is then 
that there is little evidence for black holes between 20 and 10^ Mq (Barth etal.'2004' 2005 , 2008a b , 2009 ; Greene et al. 



Greene & Ho.2007ajb 



Greene et al. 2008 ), and it appears, that there are many black holes between 10" Mq and 10' Mq, but 



probably not enough to fill the distribution. These stellar black holes add a factor to the cosmological density of black holes overall. 
Here we add some more arguments to the case that filling this hole in the black hole mass distribution is fraught with challenges: 
The finding of few black holes in this gap region is consistent with many other arguments (Kording et al.|[2002 1, saying that 
there may be no or very few intermediate mass black holes. On the other hand, since the accretion statistics of any hypothetical 
intermediate mass black holes are not known, the observed distribution of candidate sources ( |Kording et al.|[2002] l might just be 
compatible with our statistics. If we were to require that this energy density stays at the level of the 10"^ fluctuations in the microwave 
background, then MBH,mm would have to be beyond the upper end of the distribution. So whatever energy gets transmitted during 
the growth of black holes to the observed distribution, it cannot have much of any effect on the low wave-number fluctuations in the 
microwave background. 

To spin this out, one could imagine that black holes grow initially fast by feeding on dark matter dMunyaneza & Biermann 2006 



|2005j , and then continue to grow mostly by mergers; in such a scheme most of the black hole feeding would have come from dark 
matter. The gravitational waves emitted by all these growth events and mergers might obey the equation of state for dark energy 
(see Marochnik et al.|2008|l. This does not look like a convincing concept, in which we have first all dark matter to be black holes 
Abramowicz et al. (see 2009[ l, and second all dark energy gravitational waves. At present this extreme speculation does not appear 
to be a viable solution. 



4. The slope of the mass function 

The shape of the mass function is a constraint on the growth of the distribution. A simple powerlaw at high masses suggests a 
self-similar process. 

(Sil k & Takahashi|1979) l have shown how to estimate the building of a mass function from a repeated merger process, building 
on yet earlier work. Here we wish to ask, whether their approach can also explain the slope which we find here. 

We assume that all black holes grow by merging with other black holes, and that baryonic accretion is just a multiplying factor, 
like perhaps always adding a factor of 1.5 in mass. It may well be a narrow distribution of fractional mass added. 

Then writing the coalescence rate as ~ {MassY they find that the resulting mass function in its power law part becomes 
~ m"^'*''^. In our context this implies that A - 2.0. This then says that the product of cross section, typical velocity, and sticking 
probability run with this combined dependence on mass. We argued in |Gergely & Bierma nn] ( [2009| l, that the cross section runs as 
(Massyi^. On the other hand, small systems like the local group, when the earlier mergers occur, have a lower typical velocity of 
order 200 km/s than larger systems like clusters of galaxies, with velocity dispersions of order 2000 km/s. This leads to a crude 
dependence estimate of (MassY^ , with a large uncertainty. The stickiness we estimate at 100 percent for the cross-section chosen, 
given a few spiral-down orbits. The combination would imply A - 1.5, compatible with what we obtain, allowing for large error 
bars. It seems possible that the merger cross section rises more steeply with mass than what we assumed here, and an exponent close 
to unity would allow a better match with data. We conclude that mergers between black holes might be able to explain the entire 
mass distribution. 

However, we have to ask whether such a growth process could operate in a very similar way for nuclear star clusters, as for 
black holes. All such merger arguments may work as well for nuclear star clusters as for black holes surrounded by stars. 

In such a picture the upper end of the distribution is just the maximum that can be reached given the density of galaxies, and the 
mass of the central black holes. 
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5. The transition in mass 



Why is there a minimum mass in super-massive black holes? 

As we argued earlier, the data and much work by Greene et al. (|2008) and others show, that there are very few black holes near to 
and below 10^ Mq. There is a variety of exploratory ideas why this is so (Portegies Zwart et al. 2004; Munyaneza & Bier mann 2006| 



2005 I. It seems plausible to assume, that the transition from massive black holes to nuclear star clusters holds a clue to solving this 
question. A possible transition from a nuclear star cluster to a super-massive star has been discussed by rSanderS|(|l970|l, and again 



by [Portegies Zwart et al. ( 2004 1, with the latter team arguing for a transition to a super-massive black hole ( [Appenzeller & Fricke 
|1972b|at|Heger et al.|2003| l. However, the results of (Appenz eller & Fricke|1972b|a| l preclude any contribution from super-massive 
stars near or above 10° Mq, since such stars explode completely due to an instability given by General Relativity, leaving no black 
hole behind, we do need a mechanism that manages to give black holes right below this cutoff. Yungelson et al. ( [2008 1 have shown 
that wind mass loss effectively competes with agglomeration, and so would limit massive stars of a few hundred Mq to below 100 
Mq; this implies that it would be difficult to get past this barrier in mass. On the other hand, agglomeration is a runaway process, 
while stellar winds are a quasi-steady process, basically limited to the Eddington luminosity; so perhaps more extreme conditions 
are required to get a serious run-away in agglomeration. 

There are a number of processes, that contribute: One is the simple momentum exchange between stars, with a time scale of 
( |Chandrasekhar|1942[[S"pitzer|1962|[T987l|Binney & TremainefTOSTj l: 



* grav 



(9) 



where cr^. is the velocity dispersion of the stars in the system, assumed to be in virial equilibrium, is the mass of the stars, 
rii, is the density of the stars, and A is the Coulomb logarithm, typically of value 20. Considering nuclear star clusters it is hard to 
see why this process by itself would lead to a sudden transition at a specific mass, although a gravo-thermal catastrophe could in 
principle do this (Spitz er|1987 1; however, this process by itself would suggest, that the lower masses become a black hole, and the 
higher masses remain a star cluster, contrary to observation. 

On the other hand, the agglomeration of stars is governed by their collision time scale, which is 



1 



(10) 



where N^, is the total number of stars, is the density of stars, and 2* is the cross section of typical stars. We need only one 
star to start a runaway coalescence, and that is why also have the factor A^*. The question is whether either of these two processes 
or a combination of the two would allow for a transition at a specific mass of a nuclear star cluster such, that a short range of 
masses is pinpointed. McMillan & Portegies Zwart (2007 ) suggest, that there is a very long mass range, in which the process of 
agglomeration can give a large variety of masses, intermediate mass black holes. Therefore their conclusion is that this would not 
lead to a relatively sharp transition. 

Some have suggested (Perets & Alexander'2008'), that in the case of a binary black hole merger of equal masses a gravitational 
rocket effect could eject black holes from galactic centers, and one could so imagine, that all lower mass black holes might form, 
but no longer be in galactic centers. In such a speculation these black holes below the transition point would exist, but be invisible. 
However, in |Gergely & Biermann (2009 1 we show that this process is unlikely to be statistically relevant. 

In summary, accepting the process of agglomeration, what could modify the conclusion of previous authors, that a variety of 
masses is formed, and in contrast narrow down the mass range for the transition mass? 

There are several avenues to consider: 

First, galaxies grow by merging, starting from some minimum size: Could this minimum size of a central black hole correspond 
to the minimum size of a galaxy? That is hard to maintain, even considering, that Gilmore et al. ( |2007 1 have identified a minimum 
mass of order 5-10'' Mq, most of it in dark matter To grow a galaxy like ours, with a central black hole close to the low mass cut-off, 
would require so many merger events, that it is hard to see that much of any connection to the minimum galaxy could survive with 
a signature, except for properties that survive in all galaxies, independent of whether they contain a black hole at their center. 

Second, in a merger process we do obtain a central spike in dark matter from the merger density profile ( [Navarro et al.|199T 



Kravtsov et al. 



1998 



Moore et al. 



1998): 



Pdm ~ +x) y where x - r/rc, a scaled radial coordinate, and jsp is of order unity; 

various variants of this formula have been discussed ( [Kravtsov et al.|1998 i. This implies a central dark matter component with a 
mass enclosed with R of Mn ,i„ ~ F?^^''' . Applying this first just to stars as well implies a radial dependence of ~ llx^'p, of 
velocity dispersion of cr* ~ x'^^"''"'''''^. The combination yields to Taggi ~ jic^^'''''^'^/^, so to an arbitrarily short time scale of stellar 
agglomeration Taggi at the center of merged galaxies. This would then suggest, that all galaxies should have a central super-massive 
black hole, and not just those above a specific mass; this is again in contradiction to data. However, including the process of massive 
star formation, e.g.,(| Bartko et al.|2009) near the center of a galaxy might require a certain minimum amount of gaseous turn-over 
by star (star formation, mass ejection by winds, and explosions), and so it is conceivable that these processes define a threshold for 
run-away agglomeration. 

Third, in a speculation on the nature of the dark matter particle's, Munyaneza & Biermann ( 2006 2005 | l, following earlier work, 
have suggested that a keV Fermion would naturally identify a mass scale from a degenerate configuration, and this mass for a dark 
matter particle is consistent with a few 10^ Mq. This mass appears consistent with a) the low mass cutoff in the galaxy distribution 
found by Gilmore et al. ( 2007 1, b) the early star formation ([B iermann & Kusenko|2006[ Stasielak et al.|2007 Loewenstein et al. 
[2009| l, and c) new galaxy data interpretation ( [de Vega & SancheZ|2009^ . For some of these applications it is important to remember. 



10 



L.I. Caramete and P.L. Biermann: The mass function of nearby black hole candidates 



that a thermal keV particle as a dark matter candidate already is excluded ( |Shaposhnikov|2007] l; it would have to be substantially 
sub-thermal. 

Fourth, if supermassive stars at all such masses do form, evolve, and explode, it is conceivable, that they do not form a black 
hole, but a naked singularity (Joshi 2007 ), which then due to quantum loop gravity effects disperses, leaving no remnant behind. 
This would have to happen just in the mass range between about 20 and ~ 10^ Mq, but neither above or below. 

Fifth, in another concept, collapsing stars might undergo a phase transition to a SUSY-state ( |Clavelli|2006| l before forming a 
black hole, which then also blows up all the star, leaving no remnant behind. 

Finally, we need to remember, that in agglomeration of stars any massive star near to or above 10* Mq becomes unstable and 
explodes, and does so on very short time scales, leaving no remnant behind ( [Appenzeller & Fricke| 1972b ^. Therefore, all stars 
reaching to just below this limit will become a black hole, and perhaps with the abundances of the early universe, those very early 
massive stars near to but below instability threshold will also become black holes rather than completely explode. This instability 
defines a threshold, and makes it obvious, why the initial black holes - in the agglomeration picture - cannot be from above about 
10* Mq. If an evolving star moves up the main sequence, steadily growing in mass, it gets ever more unstable simply due to the 
increasing fraction of pressure held by radiation. When it approaches the Appenzeller & Fricke instability, it might just collapse, 
forming always a black hole just below threshold. A proper evolutionary calculation is required to test this. 

One could imagine, that such a dark matter core would help induce a runaway agglomeration process among stars as well; both 
run-away agglomeration, and degenerate dark matter cores, or even a combination of both, seem feasible. 

In conclusion, we have no definitive solution, why a relatively sharp transition is found in data; run-away agglomeration, perhaps 
helped along by other processes, seems like a good concept to explore further 



6. Sources and transport of ultra high energy particles 

Finally, we need to know how many black holes are out there, possibly contributing to the production of ultra high energy particles. 
In many models the power of black holes scales with mass (Falcke & Biermann 1995 ), so that we will focus on the higher masses. 
Our mass function then extrapolates within the maximum redshift used, 0.025, or approximately 100 Mpc, to 



candidates, so on the sky 0.6 per degree squared; however, note that the spatial distribution is highly inhomogeneous. Going 
to only those black holes above 10^ Mq, reduces this number to one per about 10 degree squared, so an area of about 1.8 degree 
radius, somewhat less than the correlation radius determined by the correlation analysis of Auger ( [The Pierre Auger Collaboration! 
|2007| [2008 ). To reduce the number of black holes just using the criterion of mass to a sample size of a few hu ndred as used by 



Auge r, implies a mass threshold in the distribution of about 3.5 ■ 10 Mq, firmly eliminating the radio galaxy Cen A (Cappellari et al. 
|2009 t, since it has only (5.5 ± 1.0) ■ 10^ Mq. A second criterion is required to provide a really narrow selection, such as strong radio 
emission could be (see [ Biermann et al.|2009 1, as indicative of clear particle acceleration. 

However, since most black holes probably are active in the spin-down limit (Blandford & Znajek"1977), we can adjust the 
numbers some more: The lifetime of activity in the spin-down limit is very long, and so the power emitted in the form of a jet, and 
radiation, is only a small fraction of the Eddington limit, of order percent or even less. This entails for protons, that only black holes 
above 10'' Mq can produce particles to 10^' eV. If we allow heavy nuclei with charge of order 10, then black holes above 10^ Mq 
are allowed, but we need then either an ongoing starburst to provide the particles for injection near PeV, or a latent such particle 
population left behind at substantial energy density to allow injection; this latter option appears very difficult, since cosmic rays 
leak on rather short time scales, especially at high energy such as PeV. Therefore of all the many black hole candidates ^ 10^ Mq 
discussed here, only those with a starburst and an AGN at the same time are good candidates; otherwise only those few ^10^ Mq 
and radio galaxy signatures are candidates to produce ultra high energy cosmic rays. 

Since galaxies, and especially early Hubble type galaxies, are highly clustered on the sky, and in the universe, any subset of 
super-massive black holes presumed to be candidates of ultra high energy cosmic rays will probably be even more clustered; this 
spatial distribution might show up in the arrival direction statistics. 



7. Conclusions 

We have derived the mass function of black hole candidates in the centers of galaxies down to 3 - 10* Mq. At high mass it is consistent 
with a straight steep power-law, of slope -2 in the integral mass function. In fact, since all and everyone of the black holes, which 
we could check on, has been declared a black hole by more sophisticated measurements, we feel fairly confident, that almost all 
the black hole candidates are in fact black holes; the most likely problem with such an identification is below the considered mass 
range. 

First, we conclude that this distribution function cannot possibly continue through the gap between about 30 and 10* Mq, 
consistent with Greene et al.| (|2006); Mac carone et al.| (j2007 ); Maccarone & Servillat ( 2008| l. Coming from the higher masses it 



flattens near 10' Mq, giving a contribution to the cosmological density of only 
Qbh = 2- 10-*^°'*° (12) 

consistent with many other earher estimates. Like in some other analyses our distribution begins to flatten already below 10^ Mq 
( |Lauer et al.|2007) [Benson et al.|2007| l. 
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As a corollary we note, that |Frampton | ( [20091 ) has pointed out that black holes constitute a huge reservoir of entropy: Integrating 
the mass function to obtain the total entropy we find here shows that the entropy is dominated by the biggest black holes, and is 
quite large. However, allowing for the errors in the exponent, the entropy might be logarithmically evenly distributed among all 
masses. Therefore the entropy is essentially unaffected by the speculation, that the black hole mass function may extend to lower 
masses, and in fact, it is curious that a logarithmically even distribution in entropy seems probable to be realized in Nature. However, 
this implies that the summed entropy of all super-massive black holes is a reasonably well defined quantity, a new cosmological 
parameter, that requires an understanding. 

As a second corollary we note that the early energy input of gravitational waves from the first merging black holes approaches 
oddly the energy density of dark energy at a very high redshift, of order 50, possible since very massive stars may have formed at 
even earlier redshift, such as 80 ( jBiermann & Kusenko|2006l l. 

Second, the energy input from the growth of all these black holes is at the fraction of a percent level of the microwave back- 
ground, which may mostly be lost in expansion, but some traces should remain. If some of this energy is transferred to magnetic 
fields, beyond what is given to intergalactic magnetic field from the large scale gravitational potential wells, the effect on the scat- 
tering of ultra high energy particles could be strong. In fact, this channel might present a problem. The fraction of black hole growth 
energy going into this channel is constrained to be less than a few percent. 

Third and last, just going to a distance of 100 Mpc, or redshift of 0.025, and going to 10^ Mq our results give 2.4- lO'* black holes, 
or on the sky 0.6 per every square degree. There is some evidence that every super-massive black hole is active at an observable 
level in radio emission ( Nagar et al.|2005 ). Radio emission implies non-thermal particles, and so the speculation, that each and 
everyone might contribute also to ultra high energy cosmic rays, is not immediately discountable. More specific models, with 
testable predictions, are required to investigate this possibility. For heavy nuclei as ultra high energy cosmic rays the black holes 
would require evidence of both a star-burst and nuclear activity, such as in Cen A; for protons only very massive black holes are 
allowed, such as M87. Evidence of a particle accelerator in the form of a radio galaxy limits the choice severely. 

Allowing for all such possibilities, a priori, there are many candidates every square degree on the sky; the distribution of these 
black hole candidates is highly clustered, and this clustering may survive the magnetic scattering by intergalactic magnetic fields to 
show up in the studies of the arrival directions of ultra high energy cosmic ray particles ( |Das et al.|2008t|Ryu et al.|2008| l. 
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